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We aim at deriving an equation of motion for specific sums of momentum mode occupation num- 
bers from models for electrons in periodic lattices experiencing elastic scattering, electron-phonon 
scattering or electron-electron scattering. These sums correspond to "grains" in momentum space. 
This equation of motion is supposed to involve only a moderate number of dynamical variables 
and/or exhibit a sufficiently simple structure such that neither its construction nor its analyza- 
tion/solution requires substantial numerical efi'ort. To this we end compute, by means of a projec- 
tion operator technique, a linear(ized) collision term which determines the dynamics of the above 
grain-sums. This collision term results as non-singular, finite dimensional rate matrix and may thus 
be inverted regardless of any symmetry of the underlying model. This facilitates calculations of, 
e.g., transport coefficients, as we demonstrate for a 3-dim. Anderson model featuring weak disorder. 

PACS numbers: 05.30.-d, 05.70.Ln, 72.10.-d, 05.60.Gg 



I. INTRODUCTION 

The mapping of occupation number dynamics of non- 
ideal quantum gases (as resulting from the Schroedinger 
equation) onto Bohzmann equations is a classical topic 
in the hterature As weh- known the correspond- 

ing approaches routinely rely on approximations and/or 
conceptual assumptions of some sort (for an overview see 
App. Furthermore, they often yield rather formal 
expressions for the collision terms involving, e.g, singular 
scattering rates and a diverging number of occupation 
numbers for a description based on a discrete reciprocal 
lattice. If the singularities are regularized and the de- 
scription is transferred to a continuous momentum space, 
the collision term results as a (possibly non-linear) map 
of a scalar function onto another scalar function. For 
many purposes such as calculation of transport coeffi- 
cients, electronic lifetimes etc. eigenfunctions and eigen- 
values of this map are needed. Those are, however, only 
easy to obtain for full spherical symmetry of the model, 
i.e., isotropic scattering and fully isotropic dispersion re- 
lations. If strict spherical symmetry is absent, many ap- 
proaches resort to an approximate "transport lifetime" 
which essentially amounts to an educated guess of the 
relevant eigenvalues [1]. Presumably due to this sub- 
tleties detailed calculations of, e.g., transport coefficients 
for models featuring (massively) anisotropic dispersion 
relations are rare in the literature. 

In this paper we thus suggest an alternative construc- 
tion of collision terms from pertinent underlying quan- 
tum models. These collision terms directly result as fi- 
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nite dimensional, non-singular rate matrices. Their con- 
struction as well as their numerical inversion may be per- 
formed on standard desktop computers. No singularities 
arise in the derivation of these collision terms. The ap- 
proach is based on a (optimizeable) "coarse graining" in 
momentum space. The above rates may be interpreted 
as electronic transition rates between these grains. A 
carefully chosen graining may yield a particularly favor- 
able form of the collision term, regardless of any spher- 
ical symmetry, etc., as will be demonstrated below by 
the computation of diffusion coefficients for an Ander- 
son model, cf. Sec. |Vl Our approach is an implemen- 
tation of the time-convolutionless projection operator 
method (TCL). The latter is routinely used to find an 
autonomous equation of motion for a set of variables of 
interest. Those are here essentially the accumulated oc- 
cupation numbers of the grains. Within the frame of 
the TCL approach no possibly ill-controlled, additional 
assumptions (molecular chaos, full factorizability, etc.) 
have to be employed. However, in the paper at hand we 
concentrate on a leading order truncation of the TCL- 
expansion (cf. ([1])). (Whether or not this will yield suf- 
ficiently accurate results is a somewhat subtle question 
but can in principle be checked by explicit evaluation of 
higher order terms [5|.) For a discussion of the relation 
of the present results to the afore mentioned traditional 
approaches, see paragraph below p^ . 

Our paper is organized as follows: First (Sec. |TT| we 
give a very brief introduction to the TCL projection oper- 
ator formalism. In Sec. IIIII we introduce some archetypi- 
cal condensed matter models which are routinely used to 
investigate non-ideal electronic quantum gases as arising 
from elastic scattering, electron-phonon interaction and 
electron-electron interaction. We furthermore specify our 
projection operator which is especially constructed to in- 
vestigate the dynamics of the above mentioned, occupa- 
tion number related variables. Then (Sec. IIVI) we per- 
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form leading order calculations on these models. Those 
result in general expressions for the transition rates which 
may be further evaluated on the basis of concrete disper- 
sion relations and scattering potentials. This is done in 
(Sec. |V| for a weak disorder 3-d Anderson model. It is 
thereby demonstrated how the graining can be chosen to 
cast the collision term in a form for which a relaxation- 
time approximation actually becomes exact. This facili- 
tates the simple calculation of an energy-dependent dif- 
fusion coefficient. The latter significantly differs from a 
result obtained from an approach based on full spherical 
symmetry. 

II. TIME-CONVOLUTIONLESS PROJECTIVE 
APPROACH TO RELEVANT DYNAMICS 

In this paragraph we give a short introduction to the 
TCL projection operator method ^-iSi]. In general, the 
latter is a perturbative projection operator technique 
which describes (the) reduced dynamics of a quantum 
system with a Hamiltonian of the type H = Hq + XV , 
where XV has to be small in some sense. It produces 
autonomous equations of motion for the variables of in- 
terest ("relevant information"). In order to apply this 
method one first has to construct a suitable projection 
(super) operator V which projects any density matrix 
p{t) onto its relevant part. "Relevant part" here implies 
that Vp(t), in spite of being significantly less complex 
than p(t), should still contain all variables of interest. 
Furthermore, V has to satisfy the characteristic trait of 
a projection operator V"^ p{t) = V p{t). For initial states 
with 'Pp{0) — p{Q) the TCL scheme yields a closed time- 
local differential equation for the dynamics of Vp, 

-vp{t) - mvp{t), m = ^ vjc^it), (1) 

i—l 

where the perturbative expansion used in the last equa- 
tions is in principle (formally) exact. In this paper we 
exclusively focus on initial states which satisfy the above 
relation. For a discussion of the legitimacy of this ap- 
proach see (8l-[l]|. 

As already mentioned in the introduction, we further- 
more focus on a description to leading order of ((T]). which 
is typically and certainly in our case the second order, i.e., 
we have to determine IC2 (t) . A widely accepted indicator 
for the validity of the truncation is a clear timescale sep- 
aration between the resulting relaxation dynamics and 
the decay of the correlation function, the latter being 
introduced below. 

In the literature @ one finds 

IC2{t) = I dhVC{t)C{h)V . (2) 
Jo 

where £{t)p = —ji[V{t),p] corresponds to the Liouville- 
(super-) operator. Here and in the following all equations 



are denoted in the interaction picture. For a concrete 
application we have to specify the underlying quantum 
model and further a suitable projection operator. 

III. STRUCTURE OF THE MODELS AND A 
PERTINENT PROJECTION OPERATOR 

The systems we discuss here may in general be all sorts 
of quantum gases, but for clarity and brevity we focus 
on the "spinless fermions"- type here. (This refers to 
the particles for which the collision term is to be con- 
structed.) 

^ = ^£k4«k+^ (3) 

k 



ajj.,ak are electronic creation/ annihilation operators in 
some momentum eigenmodes and Ek denotes the cor- 
responding dispersion relation. The latter depends on 
the underlying model which is assumed to describe the 
non-interacting electrons, e.g., a free electron gas, a tight 
binding model, etc.. An adequate dispersion relation may 
also be determined from an advanced solid state method 
such as density functional theory, etc.. 

Here, V refers to different pertinent interaction types, 
which we specify and investigate below in paragraph II VI 
In detail, these are elastic scattering, electron-phonon- 
interaction and electron-electron-interaction. All of them 
are treated as small perturbations (in the sense of the 
TCL method). 

Now, in order to introduce a pertinent projection op- 
erator, we firstly need to define some basic operators. 
For the non-interacting many-particle system we may di- 
rectly write down the wavenumber (momentum) depen- 
dent "single particle equilibrium density operator" on the 
mode j as: 

pT ^)«1«J + (1 - ' (4) 

with fj{p,T) = (exp((e(j) - p)/kBT) + being the 
Fermi distribution. Since we are interested in low tem- 
perature regimes we may substitute (approximate) the 
chemical potential p by the Fermi energy sf- Further we 
abbreviate f^{eF,T) as /j. 

The equilibrium density operator, again for the non- 
interacting case (, i.e.. Ho), of the total system, p^"^, may 
be written as the tensor product of the single particle 
density operators, i.e., 

P'^'^:=(g)pr- (5) 

i 

We should mention here for later reference that, while 
p'^'^ is strictly speaking just the equilibrium state of 
the non-interacting system, it is routinely considered to 
describe the equilibrium single particle properties of the 
weakly interacting system more or less correctly. Thus, 
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if single particle observables relax towards equilibrium 
due to the interactions (scattering), we expect them to 
relax towards values corresponding to 

Furthermore, the deviation of the mode occupation 
number nj = ajaj from its thermal equilibrium /j may 
be described by an operator Aj which we define as 

In order to concretize the projector further we must 
now introduce the afore mentioned coarse graining in mo- 
mentum space. Eventually this means that we have to de- 
fine domains in momentum space; we label those regions 
by Greek indices (k,??). Different occupation numbers 
corresponding to the same region will no longer be inves- 
tigated separately from each other, the remaining vari- 
ables will simply be the sums over occupation numbers 
belonging to the various grains. The concrete choice of 
the grains substantially influences the result of the follow- 
ing considerations, as will become clear in the remainder 
of this paper. However, at this point we simply math- 
ematically deflne operators A'' describing the summed 
deviations (from equilibrium) of all occupation numbers 
belonging to a common grain 

A'^^^Aj. (7) 

Furthermore, we deflne operators 

^^^= TrCfin =^^'>^^-"^"l) 

corresponds to a diagonal product operator that 
equals the equilibrium density operator on all modes ex- 
cept j, and essentially the occupation number operator 
(up to some substracted unity operator) on the mode j. 
The above deflnition of ensures that the important 
relation 

Tr{DjAi} - (9) 

is fulfilled, which turns out to be crucial for the idempo- 
tency property of the below defined projection operator. 

Moreover, we abbreviate the time dependent expecta- 
tion value of A*^, which we are mainly interested in by 
d'^{t) :=Tr{A''p(i)}, where p{t) is the density operator 
which describes the actual state of the system. 

With these preliminary definitions we construct a suit- 
able projector as follows: 

Vpit) = Tr{pit)}p^^ + E 1^(E ^J)^"(*) (10) 

Here A^^ denotes the number of momentum states in 
grain k. In the following we exclusively consider states 
satisfying Tr{/9(i)} — 1. (Of course this trace is invariant 
under unitary time evolution.) 



In general, the relevant part of the density matrix 
V p(t) is a non-factorizing state. Due to relation ([SJ the 
addition of a term containing Z?j to the equilibrium state 
p°°' in the construction (jlOp may be interpreted as adding 
a linear deviation to the equilibrium occupation number 
of mode j, while leaving the other occupation numbers 
unchanged. 

Note that a similar projector has recently been used 
in the context of investigations on electronic lifetimes in 
aluminum [T2| . In that case one is interested in describing 
the decay of a single electronic excitation and therefore 
the corresponding projector is constructed to project only 
onto one single occupation number as dynamical variable, 
whereas the projector used in this paper keeps all 
grain occupation numbers as time dependent variables 
thus giving a more detailed picture of the dynamics. In 
so far the considerations in this paper extend the previous 
work from T^]. Additionally, the are here defined in 
a slightly different way which leads to a simplification of 
the upcoming calculations. 

By exploiting Q one may straightforwardly prove that 
this projector indeed features the crucial idempotency 
property of a projection operator, i.e., = V. It fur- 
thermore obviously captures the dynamical variables of 
interest, namely the . 

Before we eventually concretely apply ([1} we make the 
following approximation for an expression that appears 
in the computation of ([T]): 

'-[V{t),p^^]^Q (11) 

The neglected commutator term essentially describes the 
dynamics of the equilibrium state of the non-interacting 
system. Eventually we are interested in a single particle 
observable. As already mentioned above, the equilibrium 
state of the non-interacting system is believed to reason- 
ably describe single particle observables in equilibrium 
even for weakly interacting systems. Since an equilib- 
rium state is constant, the above commutator should 
not significantly contribute to the relevant dynamics, 
thus we drop it. Keeping the term and performing all 
following steps eventually yields an expression which can 
be explicitly shown to be indeed negligible in the weak 
coupling limit. For clarity and briefness we omit this 
calculation here. 

Now, we explicitly evaluate ([TJ to leading, i.e., second 
order using the above projector. In order to extract from 
this equation of motion for operators an equation of mo- 
tion for the scalar observables of interest (which are the 
d'^{t)) we multiply by A'' and take a trace: 

^Tr{A^rp{t)} - Tv{A^IC2{t)Vpit)} (12) 

Exploiting dH) and [V{t),Ai^] = [V{t),n]^], furthermore 
employing dTTT) and some invariance properties of the 
trace, we find after a lengthy but straightforward cal- 
culation: 
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-TY{Di[V{h),[V{t),nu]]}d''{t) 



iGK, 
ker; 



(13) 

This is our first main result. Obviously ([13]) may be 
interpreted as a rate equation for the dynamics of the 
coarse grained occupation numbers. It thus corresponds 
to a linear(ized) collision term. The rates, which are 
typically finite and directly computable, are given by 
Rrjuit) '■— Jq dtiCrji^Xt, ti). They may in general be time- 
dependent which does not fit into the standard picture. 
But, as they are given by integrals over correlation func- 
tions, they can be expected to converge to constant val- 
ues after some correlation times under some rather mild 
conditions on the model and the graining (, cf. also para- 
graph below (fT5)) ). For a concrete example for the eval- 
uation of the rates see SecfVl 

For most routinely considered interactions (and all in- 
teractions analyzed in this paper) the particle number 
within the specifically addressed system (e.g., electron 
system) {N = J2k "k) is conserved. Thus 



ni 



(14) 



ke?7 



holds. Exploiting this feature, we may determine the 
diagonal rate terms of the rate equation from the 
non-diagonal terms by inserting (ITil) into (IT^ 



IV. APPLICATION TO MODELS FEATURING 
DIFFERENT INTERACTIONS 

A. Elastic Scattering 

Firstly, we apply the method introduced above to a 
quantum gas model according to which, in addition to a 
periodic crystal lattice potential, the electron is subject 
to a weak, non-periodic potential. The latter of course 
induces the scattering. (One may think here for example 
of a few impurity atoms in the lattice of a metal or of 
an Anderson model with very weak on-site disorder, way 
down below the critical disorder, cf. [l^. The latter 
will be analyzed in detail in Sec. |Vl) A corresponding 
Hamiltonian is given by 



k,q 



Wiq) 



n 



'■k+q 



Ok 



(17) 



Here M^(q) simply denotes the corresponding spatial 
Fourier component of the impurity potential and fl corre- 
sponds to the total number of discrete (quasi-)momenta, 
i.e., il scales with the volume of the crystal. The above 
"interaction term V" has to be inserted into Ea. dT^ to 
determine the scattering rates Rrj^it). We evaluate the 
traces in (|16p by identifying all contributing diagonal 
terms, i.e., terms that contain equally many creation and 
annihilation operators acting on a single k-mode. In do- 
ing so it proves to be advantageous that both Di and 
factorize onto the single k-modes and, furthermore, are 
diagonal on each k-mode. After some straightforward 
calculation we finally obtain for the rates 



Rr, 



E/ ^''^^ 



(15) 



Thus, the rate equation may be classified as a master 
equation which is consistent with its interpretation as a 
collision term in a Boltzmann equation. 

For the following calculations it turns out to be conve- 
nient to reformulate the trace term from (1131): 



(16) 



^TT{Di[V{h),[V{t),nu]]} = 
-{Ti{V{t)V{h)ni,Di} + c.c.) 
-{TT{V{t)nkV{h)Di} + c.c.) 



where we have used that Di commutes with rik. Note 
that the Hermitian conjugate is of identical form, respec- 
tively, except exchanged time arguments. 

The physical implications of (jl3p may become more 
transparent if it is applied to concrete models and some 
pertinent interactions are inserted. This is done for some 
examples in the next Section. 



R7jK.it) 



t 

dr 



|M/(k-i)|^ 



■cos[i(ei 

h 



£k)T] 



(18) 

with r := t — t\. (Here we have used VF(q) — W*(— q), 
since the interaction potential should be real in configu- 
ration space). 

As multiply announced, the transition rates between 
the grains result as non-singular, real numbers which are 
defined at any time . Their concrete form obviously de- 
pends on the model as well as on the graining. Rates may 
become time- independent on a timescale cx fi/A£', Ai? 
being the energy range spanned by a grain. This together 
with the requirement of the dispersion relation being ap- 
proximately linear within the grains sets the timescale 
beyond which a description based on constant rates may 
apply. This timescale may obviously be substantially 
shorter than the "infinite time limit" which is frequently 
addressed within the framework of other approaches (see 
App. If this timescale is on the order of or shorter 
than a (hypothetical) relaxation time (which one would 
obtain from this consideration) the "weak coupling limit" 
is violated, hence no description whatsoever in terms of 
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constant rates will apply. Furthermore, only rates be- 
tween grains located within the same energy range will 
be non-vanishing beyond the above timescale. Rates de- 
scribing transitions into a certain grain tend to increase 
with the corresponding cell containing more and more 
momentum modes. All of the above appears comparable 
to transition rates as obtained from Fermi's Golden Rule. 
Indeed, in the afore mentioned "long time limit" the rates 
from assume the same values that one would get 
from, e.g., a "Fermi's Golden Rule approach" (see Sec. HI 
App. \^ after regularizing the singularities and expand- 
ing the continuous occupation number function in terms 
of orthonormal "box functions" which correspond to the 
grains. So in a sense our projective approach "integrates" 
the above steps of a traditional approach and puts the 
result on an alternative theoretical footing. This "short- 
cut" may simplify detailed calculations based on the col- 
lision term significantly, as demonstrated in Sec. [V] 

Note, that the resulting rates for this model are com- 
pletely independent of any equilibrium occupation num- 
bers, i.e., the linear dynamics do not depend on the 
position in momentum space at which particles have 
been added or taken away. Or, to rephrase, deviations 
from equilibrium relax all in the same way regardless of 
whether they occur above, below, or at the Fermi-level. 



B. Electron-Phonon-Interaction 

One important mechanism that is commonly believed 
to mainly control electronic transport in metals (at 
higher temperatures) is electron-phonon-scattering. To 
investigate this case, we routinely assume a Hamiltonian 
of the following form: 



H = ^ena1an+Y^0Jiblbi + C^ ^J^^ Qk+q"k&q + /t-C.), 



k,q 



(19) 

where bl,bi are (bosonic) creation/annihilation opera- 
tors corresponding to some phonon eigenmodes (labelled 
by i) and Wi denotes the corresponding phononic dis- 
persion relation. This Hamiltonian is of the quantum 
gas type we have introduced above (jS]), but with an ad- 
ditional phononic Hamiltonian i?o,ph- The interaction 
V may be viewed as representing processes in which an 
electron is scattered under the annihilation/creation of a 
phonon such that the complete momentum is conserved. 

We choose here to project onto deviations from equi- 
librium in the electronic system only, i.e., phononic oc- 
cupation numbers are not treated as dynamical variables 
but their equilibrium values rather enter as parameters. 
To this end we keep the projection operator essentially 
as given in ([10]) but multiply p°'^ and Dj by p'^'^' the 
latter being the equilibrium state of the phononic system 
(which in itself is an ideal gas). A corresponding calcula- 
tion (which is essentially very similar to the previous one 



for electron-defect-scattering) yields 

-Tr{Di[V{h),[V{t),n^]]} = 
|W(i-k)|2 



.1 



-(/k5i-k + (i-/k)(i + .gi-k)) 



•2 cos[-(ei - £k - ^i--k){t - h)] 
n 



\W{k~i) 



-(/k(i + 5k-i) + (i-/k)gk-i) 



■2 cos[-(ei - Ek + t^k-i)(i - ^i)] 



(20) 



where we have used Tr{p<=i' ^"^bl^b^} = g^. To the corre- 
sponding rates a statement which is essentially analog to 
the one below ([T8|) applies. Again, in the long time limit 
this rate matrix assumes a form comparable to the one 
which one would obtain from a traditional approach us- 
ing the above mentioned "grain-corresponding basis func- 
tions" to represent the collision operator. However, in 
this case this only holds if the traditional approach is 
based on a linearized collision operator. So again the 
present approach integrates the steps of the traditional 
approach, including linearization. 

In many pertinent systems, the phonon energies are 
assumed to be small compared to the electron energies. 
Thus, neglecting the phonon energies (wk-i ~ 0) , we 
may further simplify (j20|) and obtain for the transition 
rates of Eq.(in]) (k ^ -q) 



-Ri)k(^) 



E 

iGK, 

ke») 



(l+5k- 



-gi_k)cos[-(ei-£k)T]. 
n 



(21) 



Thus, according to this approximation, the rate matrix 
is essentially of the same form as in the case of impurity 
scattering, except for the equilibrium phonon occupation 
numbers. Such a formulation provides (in contrast to 
other approaches) a directly and easily accessible starting 
point for massively non-isotropic systems for which any 
approximation based on full spherical symmetry must 
necessarily fail. The calculation of electronic diffusion co- 
efficients in atomic wires as performed in [14] may serve 
as an example. (For another time-convolutionless master 
equation approach to the dynamics of coupled electron- 
phonon systems see [TBj.) 



C. Electron-Electron-Interaction 

In this paragraph we apply the above method to a 
system featuring electron-electron-interaction. A corre- 
sponding Hamiltonian may read 



H 



^ . 1 ^ iy(q) 

ka,q 



Lqflkai . (22) 
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The evaluation of the trace term (1161) is basically similar 
to the previous examples. Nevertheless, identifying all 
contributing terms proves in this case to be somewhat 
subtle. However, simply following the scheme, we arrive 
after some lengthy calculation at 



(Re(VF(i-k))-Re(W^(k-l)))2.ij 
(/i(l - /i-k+i)(l - ./k) + /k/i-k+i(l - /i)) 
cos[-(ei + ei - £k - £i-k+i)T] 

~(Re(W^(l-i))-Re(T^(k-l)))2.ij 
(/k(l - /i+k-i)(l - /i) + /i.A+k-i(l - /j)) 
cos[-(ek + £i - ei+k-i - £i)t] . 



(23) 



To those rates a statement analog to the ones below 
(dHI and (I^Dl) applies. The present result is comparable 
to (and in a sense integrates) all steps of a traditional 
approach including linearization. The explicit possibility 
to choose a coarse graining in momentum space may also 
help here to cast the collision operator into an especially 
suitable form. We intend to use p3| in order to deter- 
mine transport coefhcients of one-dimensional models of 
interacting fermions as considered in [l6l | in a forthcom- 
ing paper. While in more realistic models the relevance of 
electron-electron scattering for transport has to be ques- 
tioned, its primary relevance for lifetimes of excited "hot 
electrons" in metals is undisputed. An approximation 
for the decay rate of such a hot electron results if the 
transition rates from one grain to all other grains in (pn)) 
are summed over. This turns out to be equivalent to 
employing a projector which keeps only the occupation 
number deviation within one grain as a relevant observ- 
able. Such an investigation has been performed (on the 
basis of a projector that structurally slightly differs from 
the one used in the present paper) in (l2l|. There reason- 
able lifetimes for hot electrons in aluminum have been 
found. 



Hamiltonian has the form of Eg. dTTj) . The corresponding 
dispersion relation is given by 



Ek = — 2 J(cos kx -f cos ky -I- cos k^) 



(24) 



where J is an energy determining the bandwidth 
(bandwidth=12 J). Thus the underlying ideal quantum 
gas corresponds to electrons in a simple cubic lattice. We 
focus here on random uncorrelated on-site disorder, i.e., 
the W^(q)'s are assumed to be given by independent ran- 
dom numbers, generated according to some distribution 
(box distribution, gaussian, etc.). As explained below, 
for this setting only the mean square of the l/l^(q)'s even- 
tually enters the rate matrix, thus for our purposes we 
do not even need to explicitly name the concrete type of 
distribution here. 

At this point, we specify the graining in detail. The 
momentum space is firstly divided into domains corre- 
sponding to energy shells of equal energy width /S.E (ac- 
cording to the dispersion relation). The alignment of the 
graining with surfaces of constant energy insures that 
only rates between grains from the same energy shell are 
non- vanishing, cf. paragraph below (|18|) . Subsequently 
we further divide the energy shells (labelled by E) along 
the energy gradient into grains, such that each of 
these grains contains equally many states N which im- 
plies equal volumes in momentum space. Fig. [T] shows 
an example of this type of coarse graining for the 2- 
dimensional case. We assume that the grains are small 
enough that the dispersion relation may be linearized on 
each grain. Further, we define 



1^2 



1 



(25) 



Due to the disorder being uncorrelated we get W^^ = 

— const., i.e., VF^^ is independent of the grains. Ex- 
ploiting this (and as routinely done the properties of the 
sinc-function) we can approximately perform the inte- 
gration in (jlSp . thus finding for times larger than the 
correlation time (long time limit) the non-diagonal ele- 
ments of the rate matrix (77 ^ k) 



V. APPLICATION: DIFFUSION COEFFICIENT 
OF THE 3-D ANDERSON MODEL FEATURING 
WEAK DISORDER 

In this section we aim at finding the diffusion coef- 
ficient for a 3-d Anderson model featuring very weak, 
uncorrelated disorder, i.e., such that almost all states are 
non-localized, cf. Sec. IIV Al To this end we first deter- 
mine a pertinent scattering rate matrix. Especially, we 
choose a specific coarse graining that yields a rate matrix 
which strictly obeys a certain type of relaxation time ap- 
proximation (cf. [17]) and therefore allows for a simple 
computation of a diffusion coefficient. The underlying 



^n>^ - -^—-^^E{ri),E{K) , (26) 

where Se{ij)^e{k,) is 1 if 77 and n belong to the same en- 
ergy shell and otherwise. So, as already mentioned 
only grains from the same energy shell are coupled by 
transition rates. Furthermore, due to the specific coarse 
graining, these rates are equal for all transitions within 
one energy shell, i.e., there is no specific dependence on 
r] and k. 

The total number of states within some energy shell 
E, which we label by Ne, is NgE = Ne- As explained 
above, the diagonal terms of the rate matrix may be ob- 
tained from the master equation property fEg. dlSI) ). By 
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FIG. 1: Sketch of our specific coarse graining in momentum 
space for a 2-d Anderson model. Contour lines of equal en- 
ergy define energy shells of equal energy width A_E. Within 
two exemplary shells the full detailed graining is displayed: 
the shell is orthogonally partitioned into qe domains of equal 
area. This amounts to an equal number of states per domain. 



means of the abbreviation 



h AE 



the rates may be eventually written as 



gE 



(27) 



(28) 



To repeat, the non-diagonal elements of this rate ma- 
trix are all equal and thus, this specific graining yields 
symmetric detailed balance. As well-known the equi- 
librium state, i.e., equilibrium set of "grain occupation 
probabilities" (mathematically: the eigenvector belong- 
ing to the eigenvalue 0) for this type of master equation 
is the uniform distribution within one energy shell. So 
the "normalized equilibrium vector of probabilities" cor- 
responding to some energy shell E reads: 



pE . pE _ 
^0 ■ ^0,K ~ 



1 

9e' 



(29) 



for K belonging to the energy shell E. Its elements cor- 
responding to all other shells are zero, of course. Fur- 
thermore, for this special type of matrix, all other vec- 
tors of occupation probabilities X, which are orthogonal 
to equilibrium vectors {XP^ = 0) but fall entirely into 
the respective energy shells E, span the (highly degen- 
erate) eigenspaces of the rate matrix with eigenvalues 
—R^. Thus all deviations from equilibrium at energy E 
relax exponentially with rates R^ . This is an accurate 
implementation of the scenario which is assumed to hold 
approximately if an "energy dependent relaxation time 
approximation" is performed 



depe 



Having found this rate matrix as well as its eigenvec- 
tors and eigenvalues we are all set for the calculation 
of the corresponding diffusion coefficient. There are sev- 
eral approaches to the derivation of a diffusion coefficient 
from a linear (ized) Boltzmann equation. Here we follow 
[Tsj and references therein by using the formula 



De = -u,,i?,,^w„Fo,K 



(30) 



where -D_e is the diffusion coefficient at energy E, w,;,Wk 
are the x-components (if diffusion in x-direction is con- 
sidered) of the velocities corresponding to the respective 
grains and is the (pseudo-)inverted rate matrix (ne- 
glecting the null-space). For the velocity components Uk 
we routinely plug in the slope of the dispersion relation 
in, say, x-direction at the respective grains k (group ve- 
locity). These slopes are approximately constant within 
one grain, given one employs the above mentioned grain- 
ing. Due to the symmetry of the dispersion relation the 
vector Vk,P^i^ does not have any overlap with the null- 
space of the rate matrix, i.e., Pok.'^^k.Po.k ~ 0- Hence 
R^^ from (|5D|) may simply be replaced by (— l/i?-^)(5^K, 
i.e., the inverse of the eigenvalues of the eigenspace that 
are complementary to the null-space on the correspond- 
ing shells. Hence, we may evaluate formula ([5(7| and 
eventually find for the diffusion coefficient at energy E 



k<£E 



(31) 



Thus the (shell-specific) diffusion coefficient is completely 
determined by the total decay rates R^ on the_shells and 
the averaged squared velocity in x-direction . This ex- 
pression pip may easily be evaluated using any standard 
computer. The result is displayed in Fig 15] 




FIG. 2: Energy dependent diffusion coefficient for a 3-d An- 
derson model with weak disorder according to formula pip 
(dotted curve). The structure is much richer than the one 
obtained from a free-electron approximation (solid line). 

Obviously the spectrum of diffusion coefficients fea- 
tures a non-trivial structure. Note that the maxima of 
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the diffusion coefficients are not located where the high- 
est x-components of the velocities appear, namely in the 
middle of the spectrum, but more at the edges. In this 
model the scattering potential features full spherical sym- 
metry, the only anisotropy stems from the dispersion rela- 
tion. Nevertheless an approximation based on full spher- 
ical symmetry yields unacceptable diffusion coefficients 
for almost the entire spectrum: A "free-electron approx- 
imation", i.e., replacing the true dispersion relation by a 
pertinent parabolic one (employing some effective mass), 
makes the calculation even easier (can and has been done 
analytically [l9|), but describes the diffusion coefficients 
based on the true dispersion relation only at the outer 
edges of the peaks more or less correctly, cf. solid line 
in Figl2l This insufficiency of spherical approximations is 
wellknown in the literature and various elaborate ways to 
cope with it have been developed (20l - |22j . With this An- 
derson model example we intend to demonstrate that the 
concepts developed in this paper add another straightfor- 
ward approach to anisotropy. 



dimensional systems controlled by particle-particle inter- 
actions may be addressed (cf. |16i])- If eventually conduc- 
tivity rather than the diffusion coefficient is the quantity 
of interest, diffusion coefHcients may be converted to con- 
ductivities using a generalized Einstein relation (cf. 23] ) . 
Thus the results presented here may be used as a basis 
for a range of concrete transport-theoretical investiga- 
tions. The feasibility of such an approach is supported by 
the fact that a similar approach has been demonstrated 
to yield reasonable results for electronic lifetimes in alu- 
minum [12,] . 
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VI. SUMMARY AND OUTLOOK 



Appendix A: Traditional Approaches to Quantum 
Boltzmann Equations 



We computed linear(ized) master equations describing 
the dynamics of coarse grained occupation numbers for 
electrons in periodic lattices experiencing i) elastic scat- 
tering, ii) electron-phonon scattering and iii) electron- 
electron scattering. Our approach is based on the time- 
convolutionless projection operator method. The result- 
ing rate matrices are of finite dimension and involve only 
finite rates. Their derivation does not require the regular- 
ization of any singular rates. The concrete form of those 
rate matrices, including their dimensions, depends cru- 
cially on the way the coarse graining is performed (and 
on the underlying model, of course). They may in gen- 
eral efficiently be numerically evaluated using moderate 
computing power. A suitable graining may cast the rate 
matrix in a form that facilitates further considerations 
which are based on the inversion of the rate matrix such 
as, e.g., computation of diffusion coefficients. Those in- 
versions may be computed also with moderate comput- 
ing power, regardless of any symmetry of the underlying 
model. As an instructive example we presented the cal- 
culation of the diffusion coefficient for a 3-d Anderson 
model with weak, uncorrelated disorder, finding severe 
deviations from results based on spherical approxima- 
tions. This technique to find convenient representations 
of collision terms may also serve as a basis for compu- 
tations of diffusion coefficients in other systems. One 
may think of electronic diffusion in bulk metals. Perti- 
nent (anisotropic) dispersion relations and scattering po- 
tentials could be provided by well-developed numerical 
methods, such as DFT, etc. and easily incorporated into 
the approach at hand. But also electronic diffusion in low 
dimensional systems such as nanowires embedded in bulk 
insulators may be addressed as no spherical symmetry is 
required 14]. Furthermore, diffusion coefficients in low 



Just to give some background information we intend 
to very briefly (and incompletely) review here the various 
suggestions for mappings of quantum dynamics of non- 
ideal quantum gases onto Boltzmann equations that may 
be found in the literature. To those ends we classify them 
as three main groups of approaches: i) "Fermi's Golden 
Rule approaches", ii) " Green's- functions approaches" 
and iii) "factorization approaches". Of course there are 
many interconnections between them but for clarity we 
keep to this scheme here. 

i) ; In this rather phenomenological approach the abso- 
lute square of the interaction matrix element that " con- 
nects" two occupation number eigenstates is boldly taken 
as the weight of a singular, classical transition rate be- 
tween those two states, which already implies a kind of 
random phase approximation 2J] . Furthermore, it is as- 
sumed that the probabilities of the occupation number 
eigenstates are and remain such that occupation num- 
bers factorize. This leads to an autonomous equation of 
motion in terms of mean occupation numbers (see also 

iplil). 

ii) : Here the issue is approached by setting up the hi- 
erarchy of equations of motion for the Green's functions. 
The higher order Green's functions are then re-expressed 
by two-point-Green's functions using diagrammatic tech- 
niques and corresponding approximations. This results 
in an autonomous equation of motion in terms of one- 
mode-two-time-Green's functions. The latter are directly 
interpretable as mean occupation numbers for the case of 
full equilibrium only. However, to leading order this ap- 
proach yields an outcome comparable to the result from 
i). For details see, e.g., [H^. 

iii) : The starting point are the Heisenberg equations 
of motion for the occupation numbers. Iterating those 



9 



w.r.t. time and assuming full factorizability and diag- 
onality of the occupation number operators at all times 
yields an autonomous set of equations for the mean occu- 
pation numbers. Again, this approach yields an outcome 
comparable to the result from i). For details see, e.g., 

0, [la 113. 

All of the above considerations are based on discrete 
momentum lattices. Within the frame of such a descrip- 
tion the transition rates become singular if the infinite 



time limit is considered, which is what is usually done. 
These singularities have to be regularized "by hand" be- 
fore the description may be transferred to a continuous 
momentum space. 

(Of course there are much more investigations of Boltz- 
mann equations for quantum systems, some explore yet 
different approaches to non-ideal quantum gases such as, 
e.g.,'31-33j, some treat other systems, e.g., [ISt i34 ,35i]. 
etc..) 
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